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INTRODUCTION 

In this note we briefly review the recent activity on a new class of cosmological models based on the string field theory 
(SFT) (for details see reviews [1]) and the p-adic string theory [2] focusing on their applications to the present Dark 
Energy phenomenon. The reader is referred to [3, 4, 5, 6] and refs. therein for more detailed analysis on the subject. It 
is known that the SFT and the p-adic string theory are UV-complete ones. Thus one can expect that resulting (effective) 
models should be free of pathologies. Furthermore, models originating from the SFT exhibit one general non-standard 
■ property, namely, they have terms with infinitely many derivatives, i.e. non-local terms. The higher derivative terms 

usually produce the well known Ostrogradski instability [7]. However the Ostrogradski result is related to higher than 
two but finite number of derivatives. In the case of infinitely many derivatives it is possible that instabilities do not 
appear. 

Contemporary cosmological observational data [8, 9] strongly support that the present Universe exhibits an accel- 
' erated expansion providing thereby an evidence for a dominating DE component [10]. Recent results of WMAP [9] 

together with the data on la supernovae give the following bounds for the DE state parameter wde = — 1-02+q Jg. Note 



that the present cosmological observations do not exclude an evolving DE state parameter w. Non-local models of the 
type one has from the SFT may have effective phantom behavior and are of interest for the present cosmology. To 
construct a stable model with w < — 1 one should construct the effective theory with the Null Energy Condition (NEC) 
violation from the fundamental theory, which is stable and admits quantization. This is a hint towards the SFT inspired 
■ cosmological models. 

In this note we shall demonstrate how an effective phantom behavior may appear, shall consider an indicative 
example along with cosmological perturbations arguing that the phantom phase in these models does not produce 
growing perturbations and shall outline a possible signature of such non-local models. 

THE MODEL 

We are interested in non-local models arising from SFT in the far asymptotic regime for large times. The starting point 
is the Lagrangian 



trxy/^gl + -t — d u Td^T + T 2 v(T) -A' . (1) 
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where we manifestly indicate four dimensions despite the stringy origin. We work in 1 + 3 dimensions with the 
signature (—,+,+,+), the coordinates are denoted by Greek indices fi,v,... running from to 3. Spatial indexes 
are a,b, . . . and they run from 1 to 3. Such four-dimensional action is motivated by the string field theory and details 
can be found in [4]. Gn is the Newtonian constant: %%Gn = l/Mp, where Mp is the Planck mass, a' is the string 
length squared (we do not assume a 1 = 1 /Mp), g () is the open string coupling constant. g^ v is the metric tensor, R is 
the scalar curvature, A' is a constant, T = &(a'n)T with 

□=D% = ^L^^^ v a v (2) 

and Dp being a covariant derivative, T is a scalar field primarily associated with the open string tachyon. The function 

^(oc'n) may be not a polynomial manifestly producing thereby the non-locality. In SFT one has ^(oc'D) = e~T a ' n 
but we keep this function general. j3 is a parameter determined exclusively by the conformal field theory of the string. 
Here we see the very important feature: we may have non-locality in the action. We stress here that appearance of 
non-localities is general feature of the SFT based models and exactly the one we are going to explore 2 . Fields are 
dimensionless while [g ] = length. v(f ) is an open string tachyon self-interaction. It does not have a quadratic in T 
term. Factor 1 /a' in front of the tachyon potential looks unusual and can be easily removed by a rescaling of fields. 
For our purposes it is convenient keeping all the fields dimensionless. 

Redefining the tachyon field Tb = f and introducing dimensionless coordinates x^ = x^j^faJ the above action can 
be rewritten as (where we omit bars for simplicity and put hereafter a' = 1) 



(3) 



Let us emphasize that the potential of the field Ty is V = — A?^ 2 + v(7&). Given an extremum of the potential V 



one can linearize the theory around it using 7& = 7b — t. This gives V = — \x 2 + v ^ T 2 +V(Tq). Action (3) being 



2*b 

lives 

linearized around the extremum of the potential of the field Tj, can be written as 

(4) 



= / j4 ^(l6^ + i T ^- A 



\16jcG n 2g 2 
where 

&{U) = (□ + l)^- 2 (D)-m 2 . (5) 
Here m 2 = v(7b)" and A accounts From SFT one has 

So 

,^ FI {a) = {U + \) e P D -m 2 (6) 

& is in fact the inverse propagator and it is natural to expect j3 < in the SFT case corresponding to a convergent 
propagator at large momenta. In this SFT example we can learn a very important thing. Assume m = 0. Then 
J^sft = ( D + \)eP n and the non-locality does not show up at all. Indeed, the poles of the propagator do not feel 
the exponent. Other way of thinking is that the exponential factor can be hidden by a field redefinition. The situation 
is dramatically different for ra^O. The propagator 1 /J^sft has infinitely many poles. This is the manifestation of 
the non-locality. Also we stress that physics would be totally different for full function & and its truncated series 
expansion because the pole structure may get modified considerably. 

We see that the most important role is played by the spectrum of the theory which is determined by the equation 

&{J) = {J + l)y 2 (J)-m 2 = 0. (7) 

We name it characteristic equation. It can be an algebraic or a transcendental one. Thus until examples are considered 
we shall keep ^"(n) general. The only assumption is that all roots are simple. Also the analyticity of the function & 
on the complex plane is important so that one can represent & by the convergent series expansion: 

J=£/ n O"and/„eM. (8) 

n=0 



2 Appearance of higher derivatives is not unique for this theory. Non-commutative theories, for instance, also have higher derivative, but these 
non-local structures are very different. 



Reality of coefficients is required by the hermiticity of the Lagrangian. Strictly speaking even the analyticity require- 
ment can be weakened but in this case one has to be careful with poles and consider only the appropriate domains of 
the argument were the series expansion does converge. 
Equations of motion are 



G MV = ^T MV = 4;tg£ /„£ (^□'Ta v D"- 1 - / T + a v n / T^D"- 1 - / T- 
So n =\ l=o V 

-guv (g pa d p n l xd a n n - 1 - l x+n l >un n - l ty-87iG N g tl yA, (9) 

J^(D)T = (10) 

where G^ v is the Einstein tensor, 7^ v is the energy-momentum (stress) tensor and G = G^jg 2 is a dimensionless analog 
of the Newtonian constant. It is easy to check that the Bianchi identity is satisfied on-shell and for & = f\ □ + /o the 
usual energy-momentum tensor for the massive scalar field is reproduced. Note that equation (10) is an independent 
equation consistent with system (9) due to the Bianchi identity. 

Classical solutions to equations (9) and (10) were studied and analyzed in the literature. The cornerstone in this 
analysis is the fact that action (4) is fully equivalent to the following action with many free local scalar fields 

S local = /^^(j^-^L^^is^d^d^+J^-Aj (11) 

Here there are as many scalar fields as many roots the characteristic equation (7) has. 

To see the equivalence first we write down the equations of motion for the latter local action which are 
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and 

UXi=JiXi, for all j (13) 
Note that 7; is a root of characteristic equation (7). Therefore any T; solves equation (10) just because 

&(n)zi = ^(Ji)Ti = 0, for all i 

As a consequence of linearity of equation (10) its solution is a linear combination T,. Any arbitrary coefficients in 
this summation can be accounted by the integration constants in each T;. Substituting x — T; in (9) one can readily 
find that equation (9) is transformed into equation (12). To complete the proof of the equivalence statement one has to 
show that x = T, is a general solution to equation (10). This was done for a large class of functions ,^(p). 

It is important to understand at this point that roots 7, as well as coefficients &'{Ji) can be complex. This should 
not spoil the theory just because all the local scalar fields are just mathematical functions. They do not have physical 
meaning. What is important is that the original field x must be real since it represents a physical excitation. The latter is 
easy to achieve. Namely, the roots 7, are either real or complex conjugate. Complex conjugate 7, would yield complex 
conjugate (up to an overall constant factor) solutions for t,-. Thus, it is a matter of choice of the integration constants 
to make the linear combination of all T, real. Moreover, coefficients &'{Ji) can be absorbed in the scalar fields by the 
field rescaling. We do not do this just because real &'(Ji) would indicate whether the corresponding field is a phantom 
or not. Also we note that the above proven equivalence is true for arbitrary metric. 

We continue with local action (1 1). Moreover, from now on we specialize to the spatially flat Friedmann-Robertson- 
Walker (FRW) Universe. The metric is of the form 

ds 2 = -dt 2 + a 2 (t) (dx\ + dx\ + dx\) (14) 

where a(t) is the scale factor, t is the cosmic time. The Hubble parameter is as usual H = d/a and the dot hereafter 
in this paper denotes a derivative with respect to the cosmic time t. Background solutions for x are taken to be space- 
homogeneous as well. Equations of motion become much simpler. 

3H 2 = AnG^[ M {Mi) (x 2 +M t xf) + 87tG N A, H = -4wG£^(Af,)if . (15) 



and 

tt + 3Hii + JiTi = 0, for all i (16) 

It was proven in [6] that cosmological perturbations are also equivalent in the free theory with one non-local scalar 
field (4) and in the corresponding local theory (11) with many scalar field. In the next Section we shall consider the 
most important example with two complex conjugate roots 

Obviously if we have two consequent real roots J\ and J 2 then {J\ ) and (J 2 ) are real and would have different 
signs. This is one the origin of the phantom like behavior in such models. Such phantom behavior, however, is not 
distinguishable from the already known examples. Complex roots /, in turn also produce the phantom-like behavior 
but this one is very different from the previous case. Indeed, suppose we have only two complex conjugate roots J\ 
and J2 = J\ in action (11). Then one can easily pass to the real components of the fields. However, such real fields 
will be quadratically coupled. Thus, the intuition based on the signs of kinetic terms and mass terms will be broken. 
New behavior is reflected in the fact that solutions are classically stable and perturbations do not grow even though the 
phantom divide can be crossed. 



EXAMPLE OF BACKGROUND SOLUTION 

The reduction to a one-field model is possible but practically is trivial. Nothing new arises because the structure coming 
from the SFT is not felt. One pair of complex conjugate roots J\=J and J2 = J* is a new situation worth the attention. 

It is obvious that if a complex number J is a root of then J* is a root of & as well. The system of Friedmann 
equations becomes: 

3H 2 = 4nG [3?'(J) (i 2 + Jx\) + &'(J*) (f 2 2 +/*T 2 2 )] +8nG N A, 
H = -4nG[&'{J)x l 2 + &'{J*)x 2 1 ] ■ 

This is a really new configuration much less explored in the past. The distinguishing feature is the complex coefficients. 
Because of this solutions for the scalar fields also expected to be complex. However, as it was already stressed these 
fields are not physical. What is physical is their linear combination and only the latter must be real. This is easy 
to achieve adjusting the integration constants, say, X\ = t|. What is interesting (but not surprising, though) one 
cannot have non-interacting fields passing to the real components. Precisely, fields will be quadratically coupled in 
the Lagrangian. This means that the usual intuition about field properties (like signs of coefficients in front the kinetic 
term or the mass term) may not work. In addition to the latter system of equations we can use the equations of motion 
for the scalar fields which are 

T1+3//T1+/T1 =0, x 2 + 3Hx 2 +J*x 2 = (18) 

Equations (17) and (18) describe the late time evolution of the model with Lagrangian (3). This model possess a 
solution with the scalar field tending to the minimum of the potential (i.e. x — > 0) and the Hubble parameter going to 
the constant. Such solution was constructed numerically and was proven to be a solution [11]. Also the asymptotic 
form of this solution was derived in [4]. The idea is to compute a solution to (18) using the constant H = Hq and then 
compute the correction to H using (17). Then the procedure can be iterated to give higher and higher corrections. It 
was proven in [4] such iterations do converge. 

Solution to (18) with constant H = Hq is obviously 

Ti - Ti+e ^ + %\-e a -\ %2 = T2+e< 1 + X 2 -e a -' (19) 



where a± = -j 4 y— 1 ± — ^71 J ■ Consider Ti. Only one term in the solution converges when t — > °° in general (if 

both converge we select the slowest one). Let's assume it is the first one proportional to Ti+ constant. Then in order to 
pick the (slowest) convergent solution we put Ti_ =0. The similar term with T 2+ will converge in T 2 3 . To assure the 
initial field x is real we have to impose T 2+ = T* + . Here we define X\+ = To and a+ = a. 



3 Here we adopt the rule i/i* = ^fz meaning that the phase of the complex number runs in the interval [— 71, Tt) and for z = re'f the square root is 



The first correction to the constant Hubble parameter in case only decaying modes in x are present is 



// = //o + /! = // + /z (t? + t 1 * 2 ) 



(20) 



Constant ho is not independent and is related with To- We note that h is of order x 2 . The last expression is a good 
approximation for H in the asymptotic regime when h <C Hq. 
For the energy and the pressure we find 



1 

P = 2 



'&'{J)x 2 (a 2 + j)+,?'(J*)x* 2 (a* 2 + J*^ 
& (J ) x\ (a 2 - J) + {J* ) rf 2 (a* 2 - J* ) 



+g 2 A, 



and consequently for the state parameter and the speed of sound 

t2 



w ; 



-1 - 



&'(j)j%\ + &'{j*)r f %? 



2 _^'(J)ax 2 (a 2 -j)+^ / (J*)a*xf(a* 2 -J*] 
Cs ~ ,^'(J)ax 2 (a 2 + J) + ^'{J*)a*x* l 2 (a* 2 + J*) 

Further one can find the scale factor to be 



a = flo ex P Hot + 



ho ( XT 



a a* 



(21) 



(22) 



(23) 



It is manifest that the phantom divide is crossed thus periodically triggering the phantom phase of the evolution. 
The next problem is the behavior of the cosmological perturbations. The main question: do they grow as one would 
expect in the ordinary phantom model or not? 



COSMOLOGICAL PERTURBATIONS IN THE NEIGHBORHOOD OF THE 
SOLUTION WITH COMPLEX MASSES 



Configurations with a single scalar field were widely studied and those appearing in the non-local models do not 
have any distinguished properties. Roughly speaking configurations with many scalar fields were explored as well but 
we have here new models featuring complex masses and complex coefficients in front of the kinetic terms. As it was 
stressed above there is no problem with this for the physics of our models while properties of such models, in particular 
the cosmological perturbations with such scalar fields were not studied in depth. Thus we focus on perturbations in the 
configuration with complex roots J. The simplest case is one pair of complex conjugate roots where the background 
quantities were derived in previous Section. 

To find out the energy density perturbations we must use a system of equations on e which is the invariant energy 
density perturbation and on £ which is the gauge invariant scalar field perturbation (see, for instance, [12, 13, 14]). 
The system is 



(P + P ) ( C + (3ff + a + a*) C + (-3H + £ 
&'(J*)a* 2 xf - ,^'(J)a 2 x 2 ] C + 2# 2 Ae) 



r 

a* 



(24) 



(p +p)(^e + e# (8 + 3c 2 ) + e (^15// 2 + 9# 2 c; 



k 2 

2 - 2 i _„-2tfo' 



2k 2 J?'{J)&'{J*)a 2 a* 2 x 2 x^ fJ__J^ ] e2{ . Ho+a+alt , 

a a* ' * 



(25) 



where we should use 



H = 2h (r 2 a + T*V 
p + P = ,?'{J)T\a 2 + &'{J*)Tfa* 2 



Cs 2 = 



^'{J)a% 2 (a 2 -J)+&'{J*)a*xf (a* 2 -J*) 
&'(J)ax\ (a 2 + J) + ^'(J*)a*xf (a* 2 + /*) 



Here H is given by (20), a is given by (23), energy and pressure are given by (21) and w and c 2 follow from (22). 
Note that e is a real function and £ is an imaginary function of time only. The latter system of equations is ready to 
be solved numerically but in order to get some insight in what is going on it is instructive to make some assumptions 
about the value J. This makes some analytic progress possible. 

We recall the SFT origin of the model. Practically this means that values of J are determined with the string scales 
while Ho is expected to be much smaller. Therefore, it is natural to assume that \ y/J\ » Ho- This implies 

a s=s iVJ, oc 2 w J — i3H VJ 
a* w -iVJ*, a* 2 ss -J* + i3H VJ* 
under this assumption the system (24,25) becomes 

(p +p)'C + 2i (&'{J)J>clsfF- &'{J* )r t?Vj) C + 
+ (-i6ho(p+ P )(VjT 2 - Vrxf) + (p+ p )^)^ = (26) 

=-i(yj+Vj*)2giA£ 

(p + pYe+tii^'ijyVjxl-^'iJ^rVrxf)- 

■6iH V^j &\J)Jt\ + + 6iH VJ^} & (J*)J* rr 2 ) = (2?) 
where we should use 

p+p = -{^'{J)T 2 l J + &'{J*)Tfj*) 

Using the explicit expression for X\ = Toe"' S3 %o e '^' and representing a = x/2 + iy/2 one can write for any A 

At 2 +A*x* 2 = 2A e xr cos(yf + <p) 
-/(At 2 + A*Tj 2 ) = 2A e*' sin(>t + <p) 

Introducing % — i-^re^C, the equations of interest can be written as 

80^ 



*0 

= — I 

a 



cos(yf + r)% + 2(\/x 2 + y 2 sin (yt + r + <p b ) -xcos(yf + r))%+ 
cos(yf + r) (bho\/x 2 +y 2 e xt sin(yf - <p b ) + ^ +x 2 I - 



(28) 

"o / 

-2x yjx 2 +y 2 sin(yf + r + (p b f) % = -2y£ 

cos(yf + r)e + 2^Jx 2 +y 2 sin(yf + r - <p fc )e+ 

,„ /M 2 , , 2k 2 y (29) 



where all the constant coefficients are real, e and % are real, x is expected to be negative, R > 0, r depends on &'(J), 

and x y 

q>b = arcsin -^=^= , <p c = arcsin — . 

The most alarming points of the evolution are yt + r = j+nn where the coefficients in front of second derivatives 
become zero. Numeric integration may hit problems at these points if the precision is not very high. In the neighbor- 
hood of these points one has 

t%-2x + 2xx = 2e (30) 
2k 2 

te-2e-3H e = ^ r X- (31) 
a o 

For negative x the solution for e around t = is e = £o + £\t + . . . meaning that these points are not singular for the 
above system of equations. 

A typical behavior for the function e is dumped oscillations. Such a behavior does not depend on the wavenumber 
meaning that perturbations with complex conjugate scalar fields do vanish. This is different from usual models with 
real scalar fields where different regimes exist and most likely growing modes are present. 

This analysis can be easily extended to many complex roots and the main result, absence of growing modes at 
large times will remain. Indeed, the background solution is determined with the component having largest |x|. Thus 
the background quantities can be considered approximately the same. New equations (in total we have N equations if 
we have N roots) can also be accounted numerically and one can demonstrate that the main effect of absence of the 
growing perturbations in e does remain. The best way to see this, however, is to go to the exactly solvable system in 
which the analysis is more transparent. The detailed analysis of this question is to be found in the forthcoming paper 
[16]. 



COUPLING TO DILATON AND THE COSMOLOGICAL SIGNATURE OF 

NON-LOCAL MODELS 

It is natural to consider a more general coupling of the open string modes with the closed string excitations. The first 
non-trivial step is to add the dilaton to the metric. The action (1) is minimally modified as follows 

S = y ^ a ^^_|_^_ (/? + ( ^ (J)) 2_ (/((I))) + ^ ^-0/2 ^_l^ W r + If 2 ^) -^^v(f)^) -A') . (32) 

Such a modified model was analyzed in details in [5]. Here <t> is the dilaton field and U (<E>) is the dilaton potential. 
The main point is that on one hand the dilaton field is dynamical while on the other hand the restrictions on the 
rate of change of the Newtonian constant are very narrow (for example, from the Solar system measurements). The 
latter circumstance allows to use a series expansion for the dilaton field around some constant value. Effectively this 
means that equations derived in previous Sections get modified with small corrections. However, the presence of the 
dynamical dilaton field dictates more restrictions on solutions and modifies the possible vacuum expectation values 
for the field T as well. This in turn alters the spectrum of roots /, making the system significantly different from the 
the just tachyon dynamics. 

There is one exceptional case for the above system when the tachyon potential has a minimum with both first 
and second derivatives vanishing. Physically phis means that the potential is flat in the minimum. This corresponds 
to the real root J = if v(f ) ~ T 4 . In this particular case dilaton corrections render small deviations from 7 = 0. 
These corrections become comparable to the Hubble scale. It is possible because of the smallness of the dilaton rate 
of change. As a consequence all the parameters in the solution turn out to be of the present cosmological scales. 
The resulting behavior are oscillations with a period of order of 1 Gyr in the state parameter w which is detectable. 
Moreover, there are indications of oscillations of this kind in the literature [15]. 

It is the situation which one cannot easily achieve having just an open string sector while normally all the string 
quantities are of the order of a'. 



SUMMARY AND OUTLOOK 



To summarize we have put in line the appearance of non-local models from SFT, construction of asymptotic solutions 
in such models, analysis of perturbations and a possibility to extract observationally detectable effects. 

Indeed, models coming out of the UV-complete theory are interesting because they are expected to be free of 
pathologies. Appearance of non-localities is natural in SFT while their impact on cosmological models is a novel 
effect. Solutions in the non-local models in general are not constructed yet. There are numeric attempts while analytic 
solutions are known only in the linear models. 

The example of perturbations with complex roots carried out in this note demonstrates that linear perturbations 
can be confined thus not destroying the system. It is not evident from the very beginning and supports the claim 
that the SFT based models are stable. The case of complex 7, has never been studied in general and deserves deeper 
investigation. This problem and other technical questions related to the equations of the present paper will be addressed 
in the forthcoming publications. 

Another interesting question is the study of models which possess exact solutions preserving the asymptotic 
behavior. This would drastically simplify the analysis of the perturbations. Also generalization of the localization 
method to the presence of the dilaton field is a yet unexplored problem worth further investigation. 

As a more ambitious problem which is of great importance is a construction of the formalism analogous to presented 
in this paper for a model with self-interacting non-local scalar field. Such models play important role in the SFT. For 
instance, rolling tachyon dynamics is governed by action (3) with a polynomial potential of fourth degree. However, 
even background solutions are not very well understood because there is no general analytic way of solving non-local 
non-linear equations. On the other hand it follows from the present analysis that passing to a local system with many 
fields is vital for the construction of perturbation equations. 

Looking a step further it is interesting to consider perturbations in other non-local models coming from the SFT. 
For instance, models where open and closed string modes are non-minimally coupled may be of interest in cosmology. 
An example of the classical solution is presented in [5]. Furthermore it should be possible to extend the formalism 
presented in this paper to other models involving non-localities like modified gravity setups. 
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